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ABSTRACT 

We show that measurements of stellar proper motions in dwarf spheroidal galaxies provide a powerful 
probe of the nature of dark matter. Allowing for general dark matter density profiles and stellar 
velocity anisotropy profiles, we show that the log-slope of the dark matter profile at about twice the 
stellar core (King) radius can be measured to within ±0.2 when the proper motions of 200 stars are 
added to standard line-of-sight velocity dispersion data. This measurement of the log-slope provides 
a test of Cold and Warm Dark Matter theories at a sensitivity not possible with line-of-sight velocity 
dispersion measurements alone. The upcoming SIM PlanetQuest will have the sensitivity to obtain 
the required number of proper motions in Milky Way dwarf spheroidal galaxies. 
Subject headings: Cosmology: dark matter, theory-galaxies: kinematics and dynamics-Astrometry 



1. INTRODUCTION 

In a Universe dominated by collisionless, Cold Dark 
Matter (CDM), the central densities of galaxy halos are 
high, and rise to form cusps with p ~ r" 1 , as r — > 
(jDubinski fc Carlberd [i~99ll : iDiemand et alj I2005D . In 
contrast, a large class of alternatives to CDM, here gener- 
ically called Warm Dark Matter (WDM), are expected 
to have lower ce ntral densities and constant density cores 
at small radii dDalcanton fc Hoganl 120011: i Kapl inghat 
120051 : ICembranos et all 120051 : IStrigari et al.l l2006al ). In 
the majority of WDM models, phase space arguments 
predict that density cores will be more prominent in 
lower mass halos. Current observations of the small- 
est rotationally-supported galaxies are somewhat am- 
biguous. Some galaxies prefer cored halo s, some pre- 
fer cusps, and others are well fit by either (|Simon et al.1 
[20051 : iKuzio de Narav et al.ll2006h . Moreover, the effects 
of non-circular motions and other systematic issues are 
yet to be sorted out for these systems. 

The dwarf spheroidal (dSph) satellite galaxies of the 
Milky Way provide a potentially superior laboratory for 
studying the nature of dark matter. Observationally, 
their proximity (~ 100 kpc) allows kinematic studies 
of individual stars. Theoretically, their small masses 
make them ideal candidates for prominent WDM cores. 
Moreover, the CDM prediction of cuspy central density 
profiles is robust for these satellite systems because the 
cusps are stable to any tidal inter actions that may have 
occurred (Kazantz idis et all |2004|) . Unfortunately, ob- 
taining dSph density profile slopes has been difficult. 
The best current constraints on their mass profiles come 
from analyses of ~ 200 line-of-sight (L OS) velocities (e.g. 
lLokas et al.ll2005l iMunoz et al.ll2005[ ). The stellar veloc- 
ity anisotropy is a major source of degeneracy in LOS 
velocity dispersion modeling, and, as we show below, the 
log-slope of the underlying density profile remains practi- 
cally unconstrained even if the number of observed stars 
is increased by a factor of ~ 10. 

In this Letter, we examine the prospects of constraining 
the dark matter density profiles of dSphs by combining 
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stellar LOS velocities with proper motion measurements. 
We show that ~ 200 proper motions at ~ 5 km s _1 accu- 
racy can break the relevant degeneracies and determine 
the dark matter distribution in the vicinity of the dSph 
stellar core radius, r^ng- Specifically both the dark halo 
density and the local log-slope of the dark halo density 
profile at ~ 2 r^ng may be determined with > 5 times 
higher precision than from LOS velocity dispersion data 
alone. Using two well-motivated examples of a cusp and 
core halo, we show that the local log-slope measurement 
can distinguish between them at greater than the 3-<r 
level. We discuss our results in the context of NASA's 
SIM PlanetQuest (Space Interferometry Mission) which 
will have the sensitivity to measure proper motions of 
stars at this accuracy in multiple dSphs. 

Astrometry as a means to constrain the central den- 
sity profiles of dSphs was previously considered by 
IWilkinson et al.l ((2002). Using a two-parameter family 
of models for the dSphs, these authors show that proper 
motions will be enough to completely reconstruct the 
profile if the underlying shape follows their adopted form. 
However, CDM halos, and likely their WDM counter- 
parts, are described by a minimum of four unknown pa- 
rameters: a scale density, a scale radius, an asymptotic 
inner slope, and an asymptotic outer slope. By marginal- 
izing over all of these parameters as well as the velocity 
anisotropy of the stars, we demonstrate that the asymp- 
totic slopes are never well-determined, even with proper 
motions. We show that what is best constrained is the 
log-slope, 7(r*) = — d\n.p/d\n.r\ r=Tt , at a radius com- 
parable to the King core radius, specifically r* ~ 2 rking- 
We note that unlike the asymptotic r — > slope, the log- 
slope at r* ~ 0.4 kpc is known from CDM simulations. 
This radius corresponds to ~ 1% of the relevant halo 
virial radius, and is well-resolved in current simulations. 

2. MASS MODELING 

We define the three-dimensional (spatial, not pro- 
jected) components of a star's velocity to be v r , vg, and 
I),)). The velocity component along the line of sight is 
then vios = iv cos 6* + vg sin 9, where z ■ r = cos# and 
z is the line-of-sight direction. The components parallel 
and tangential to the radius vector R in the plane of the 
sky are vr = v r sin 9 + vg cos 9 and v t = v^ , respectively. 
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For each component, the velocity dispersion is defined as 
of = (v 2 ). We will assume of — of. 

The velocity dispersion for each observed component 
can be constructed by solving the Jeans equation for the 
three-dimensional stellar radial velocity dispersion pro- 
file oy(r) and integrating along the line of sight. We 
note that even in the cas e of tidally disturbed dwarfs, 
iKlimentowski et al.l (|2006h have shown that dSph veloc- 
ity dispersions are well modeled by the Jeans equation, 
as long as unbound, interloper stars are removed with 
standard procedures. We derive the three resulting ob- 
servable velocity dispersions: 
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Here (3(r) = 1 — o 2 ja 2 is the stellar velocity anisotropy, 
is the surface density of stars, and v*(r) is the 
three-dimensional density of stars. It is clear from in- 
spection that each component depends on (3 in a different 
fashion, and therefore can be used together to con strain 
its value. For h(R) and v±(r) we use a King profile (King 
Il962f) . which is characterized by a core radius, n^ng, and 
tidal radius, r t . We adopt values that describe the sur- 
face density of Draco: r t = 0.93 kpc and rki ng = 0.18 
kpc. Note that the remaining dSphs have similar King 
concentrations, r t /r^ a& ~ 5, with Sextans having the 
largest ratio r t /rki ng ~ 10. Our results do not change 
significantly as we vary r t /rking (equivalent to looking at 
different dSphs). 

In Eqs. (HJ), and ([3]), the radial stellar velocity dis- 
persion, ay, depends on the total mass distribution, and 
thus the parameters describing the dark matter density 
profile. We will consider the following general parame- 
terization of the dark matter density profile, 
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Here, the value of a sets the asymptotic inner slope, and 
different combinations of b and c set the transition to the 
asymptotic outer slope. For the sp ecific choice (a.b, c) = 
(1, 1, 3), we have an NFW profile (|Navarro et al.lll996D . 
We denote this as our cusp case below. For our core case 
we use (a, 6, c) = (0, 1.5 , 3), corresponding to a Burk- 
ert profile (Burkcrt 1995). We take these two models to 
be representative of the predictions of CDM and WDM 
models. The Burkert profile for the core case is moti- 
vated by the expectation that WDM halos will mimic 
CDM halos at large radiu s . This was seen in the WDM 
simulations of IColm et alJ (2000). The Burkert choice is 
also conservative compared to the often-used isothermal 
core with (a, b, c) = (0,2,2), which is more divergent in 
shape from an NFW and would be easier to distinguish 
observationally. Regardless, our methods are robust to 
changes in the underlying form of the density profile. 

Eq. (j4|) allows considerable flexibility in overall form, 
and the five shape parameters (a, 6, c, r$, po) are in many 
cases degenerate. However, there are a number of phys- 
ically relevant quantities that may be derived for any 



set of the five shape parameters. The first is the log- 
slope of the dark matter density profile, defined as j(r) — 
— d\n p(r)/d\nr. For the density profile in Eq. ^ this is 
given by 7(7*) = a — (a — c)(r/r ) b /[l + (r/r ) b ]. Other 
quantities of physical interest are the integrated mass 
within a given radius, M(r) 1 and the physical density at 
a given radius, p(r), which are clearly obtained for a de- 
generate set of shape parameters. Below, we show that 
while the shape parameters are not well constrained by 
dSph velocity data, the physical quantities of interest at 
the scale of the stellar core radius, r* ~ 2rki ng , may be 
constrained to high precision. 

3. FORECASTING ERRORS ON PARAMETERS 

Our goal is to estimate the accuracy with which the 
velocity components of stars in dSphs can be used to 
probe the underlying dark matter distribution. We will 
consider a model with six independent parameters: a, 
b, c, po, r s , and (3 = constant. We will consider gen- 
eralized (3{r) forms below. In order to keep the profile 
shape relatively smooth (as is expected for dark matter 
halo profiles) we restrict the range of b and c by adding 
Gaussian priors of ±2. 

The errors attainable on these parameters will de- 
pend on the covariance matrix, which we will approx- 
imate by the 6x6 Fisher information matrix F %3 — 
(d 2 \nC/d Pl d Pj ) (|Kendall fc Stuartl [19691. The inverse 
of the Fisher matrix, F~ , provides an estimate of the 

covariance between the parameters, and \J F^ 1 approx- 
imates the error in the estimate on the parameter p t . 

The Cramer-Rao inequality guarantees that y F^ 1 is 
the minimum possible variance on the ith parameter for 
an unbiased estimator. Using F _1 in place of the true 
covariance matrix involves approximating the likelihood 
function of the parameters as Gaussian near its peak, so 
F _1 will be a good approximation to the errors on param- 
eters that are well-constrained. The Fisher matrix also 
provides information about degeneracies between param- 
eters but obviously should not be trusted for estimates 
of the error along these degeneracy directions. 

We pick large radial bins to compute the velocity dis- 
persions and check that this uncorrelates the different 
bins. Then the elements of F are given by 
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The sum is over I radial bins and M refers to the three 
velocity "methods" - one line-of-sight and two compo- 
nents in the plane of the sky. The errors on the velocity 
dispersion are represented by €a«- We choose bins of 
equal width in distance, so that there are approximately 
an equal number of stars in each radial bin. As long as we 
distribute equal number of stars in each bin, the results 
we present below are insensitive to the binning scheme, 
except in the limit of very few bins, or in the limit of 
small numbers of stars per bin. 

To model the errors on the velocity dispersions, we de- 
fine e\. le — ([cr Me — (o~f ie }} 2 } ■ We assume that the errors 
on the velocity of each star are Gaussian and that the 
theory error (from the distribution function) and exper- 
imental error are summed in quadrature, 
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TABLE 1 

Fisher Matrix 1-ct errors on parameters. 



Parameter 


Fiducial Value 


Error Quoted 


200 LOS 


1000 LOS 


200 LOS 


1000 LOS 


1000 LOS 








+ SIM 


+ SIM 


+ 200 SIM 


+ 200 SIM 


+ 500 SIM 


70* = 2r Hng ) 


1.3 (0.32) 


A7* 


2.0 (2.2) 


0.91 (1.0) 


0.28 (0.24) 


0.21 (0.23) 


0.16 (0.15) 


M(r* = 2r king ) [1O 7 M ] 


1.3 (5.3) 


AMj/M t 


0.27 (0.38) 


0.12 (0.17) 


0.11 (0.11) 


0.07 (0.08) 


0.06 (0.05) 


p(r* = 2r ki ) [10 7 M Q kpc" 3 ] 


4.0 (2.5) 


Ap*/p* 


1.6 (0.73) 


0.73 (0.33) 


0.16 (0.16) 


0.13 (0.12) 


0.09 (0.09) 


Vuiax [kms -1 ] 


22 (28) 


AVmax /Vmax 


1.6 (1.7) 


0.71 (0.76) 


0.68 (0.93) 


0.45 (0.62) 


0.38 (0.52) 


P 


0(0) 


A/3 


2.1 (1.4) 


0.97 (0.62) 


0.16 (0.16) 


0.15 (0.15) 


0.10 (0.10) 


a 


1 (0) 


Ao 


4.4 (4.6) 


2.0 (2.1) 


0.78 (0.87) 


0.55 (0.66) 


0.44 (0.51) 


po [1O 7 M kpc" 3 ] 


1.0 (3.2) 


Apo/po 


20 (9.2) 


9.1 (4.2) 


5.9 (2.0) 


4.0 (1.5) 


3.3 (1.1) 


ro [kpc] 


2.0 (1.5) 


Aro/ro 


8.3 (5.1) 


3.8 (2.3) 


3.3 (1.8) 


2.2 (1.2) 


1.8 (1.0) 



Note. — Errors refer to the 1 — a range obtained from a full marginalization over the 6 input parameters (a, 6, c, ro, po, /3). Quantities 
without (with) parenthesis refer to the cusp (core) model described in the text. In the last 5 columns, LOS refers to the number of 
linc-of-sight stars used to derive the quoted errors, and SIM refers to the number proper motions. 
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Fig. 1. — The 1-a error on the log-slope of the density profile 
as a function of the radius for 1000 LOS velocities combined with 
200 proper motions. The left and right panels show our results for 
the cusp and core cases, respectively. The value of r* is set to be 
where the error is minimized. 

where n is the number of stars in each bin (taken here to 
be constant from bin to bin). Here a\ is the dispersion 
in each bin determined from Eqs. C}-©, and a™ Il rep- 
resents the measurement error in the velocities of stars 
in that bin. For LOS velocities, the spectroscopic resolu- 
tion implies errors ojJJs ~ lkms -1 , which are negligible 
compared to the underlying dispersion, c* ~ lOkms . 
SIM can achieve ~ 10/zas yr _1 proper motions for faint 
stars 3 . At 100 kpc, this translates to an accuracy of ~ 5 
km s _1 . We will take this as a typical measurement error 
for velocities in the plane of the sky. 

The derivatives in Eq. (O, and thus the errors at- 
tainable on any of the parameters, depend on the lo- 
cation in the true set of parameter space. To exam- 
ine how the errors vary as a function of the true set of 
parameters, we choose two fiducial models. Both mod- 
els produce a "typical" dSph velocity dispersion profile, 
which is roughly flat at ~ lOkrns -1 out to r t . This 
requirement does not completely fix the profile. For 
our cusp model (NFW) we add the further restriction 
that rn and pn fall within the expected CDM range 
(|Strigari et alj l2006bf) . For our core model (Burkert) 
we set tq and po to give a central phase space density 
Q ~ 10~ 5 Mp) pc~ 3 (km s~ 1 )~ 3 , which can b e produced 
in some WDM models (|Strigari et alj|2006al) . The val- 
ues of ro and po for both models are listed in Table Q] 
as are the implied log-slopes 7(7*) = 1.3 (cusp) and 0.32 
(core) at the characteristic radius r*, which is set below 
at the value 2 r k j ng = 0.36 kpc. 

3 http://planetquest.jpl.nasa.gov/SIM/sim_AstroIndex.cfm 



4. CONSTRAINING THE LOG-SLOPE 

The log-slope of the density profile varies with radius, 
and so does the minimal error attainable on it. We are 
interested in searching for the radius, r*, where the log- 
slope is best constrained. Fig. [1] shows the I — a error on 
j(r) computed as a function of radius using 1000 line-of- 
sight velocities and 200 proper motions for the cusp (left) 
and core (right) models. In both cases, the error reaches 
a minimum of A7 ~ 0.2 at ~ 2r k i ng and we explicitly 
adopt r* = 2 r k i ng for the rest of the discussion. We find 
that the value of and the minimum error depend some- 
what on the fiducial model. As shown in Fig. [U profiles 
with cusps are slightly better constrained than profiles 
with cores, and r* occurs at slightly smaller radii for 
cusped models. We also find that as r t /r k i ng increases, 
the minimum error decreases and r* increases. We have 
chosen a conservatively small ratio here, r t /r k i ng ~ 5. 

Table [1] summarizes our results. We list 1-a errors at- 
tainable on several quantities determined by marginaliz- 
ing over the 6 fit parameters discussed above. We present 
errors for five different combinations for the number of 
LOS stars and proper motion stars. In all cases, the 
errors on the "physical" parameters, M(r*), p(r+), and 
7(r*), decrease significantly when proper motions are in- 
cluded. Quantities listed without (with) parentheses cor- 
respond to the cusp (core) halo case. Note that the errors 
on the halo shape parameters b and c are essentially the 
same as the priors we set on them of ±2 and hence we 
do not list them in the table. 

A quantity that is well constrained without the addi- 
tion of proper motion information is the mass M(rv), 
which can be determined to 27% (38%) for cusped 
(cored) models with current data. The mass can be de- 
termined to a remarkable 7 — 8% accuracy when 200 SIM 
stars are added. The maximum circular velocity is rel- 
atively unconstrained by LOS data alone, but becomes 
determined to within 45% (62%) for cusped (cored) mod- 
els with 200 proper motions. More interesting for the 
nature of dark matter is the log-slope of the density pro- 
file at r*. Even with 1000 LOS velocities, the log-slope in 
the cusp case is virtually unconstrained, but with the ad- 
dition of 200 SIM stars the log-slope is well-determined, 
e.g. 7(7*) = 1.3 ± 0.2 for the case of an underlying cusp. 

The power of velocity dispersions to constrain the 
density and mass at r* ~ r k j ng is relatively sim- 
ple to understand by examining the observable veloc- 
ity dispersion components (Eqs. [IJ [2] and [3|). These 
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Fig. 2. — The 68% and 95% confidence regions for the dark 
halo density profile slope measured at r* = 2ri;i n g and velocity 
anisotropy /3. The left (right) panels correspond to the cusp (core) 
halo model and the small x's indicate the fiducial values. Upper 
panels show the errors with 1000 LOS velocities and the bottom 
panels show the errors for an additional 200 proper motions. 

observable quantities depend on the three-dimensional 
stellar velocity dispersion, which scales as a 2 (r) = 
v -i r -2i3 f rt Q Vi< (j)M{r)r 2 P~ 2 <ir oc r 2 ~ 7 * for a power- 
law stellar distribution v±(r) and constant (3. The ma- 
jority of stars reside at projected radii rking ^ R J$ ft, 
where the stellar distribution is falling rapidly v± ~ r - 3 - 5 . 
In this case, for (3 = 0, the LOS component scales as 
<rf oa {R) oc Q r -o.5-7* ( r 2 _ #2)-i/2 dr . and is dominated 
by the mass and density profile at the smallest relevant 
radii, r ~ rking- For R < r^- mg , v+ oc r" 1 and af os is sim- 
ilarly dominated by r ~ r^ing contributions. Therefore 
we expect the strongest constraints on the dark matter 
profile at ~ n<i n g. For (3^0 similar arguments hold 
but there is a manifest degeneracy between (3 and the 
log-slope 7*. As discussed above and illustrated in Table 
1, this degeneracy can be broken by including both LOS 
information (Eq. [T]) with tangential information (Eqs. 

HE]). 

With the constraints on 7*, are we able to distinguish 
cored from cusped models? To answer this question, in 
Figure [2] we present 68% and 95% error contours in the 
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two-dimensional plane of 7*-/?. This figure clearly shows 
the utility of combining both LOS and proper motions. 
Due to the degeneracy with (3, a sample of ~ 1000 line- 
of-sight stars is unable to distinguish between our two 
fiducial models, though the models can be clearly distin- 
guished with the addition of proper motions. In the /3-7* 
plane, the principal constraining power comes from com- 
bining the LOS and R velocity dispersions because the (3 
dependence in these two quantities comes with opposite 
signs (Eqs. HJ 02). 

Finally, we address the fact that (3 could in princi- 
ple vary significantly with radius. We have repeated 
our analysis assuming that (3{r) = (3 + /3ir 2 /(Y 2 + r 2 ), 
and marginalize over (3q, (3i and rp. We find very small 
changes to the 7(7*) error. For reasons similar to those 
given above, the three different velocity dispersion pro- 
files constrain /?(r = r*) to high accuracy, breaking the 
/3-7 degeneracy. 

5. SUMMARY 

We have shown that the measurement of proper mo- 
tions for ~ 200 stars in a typical dSph can be combined 
with current line-of-sight velocity measurements to con- 
strain the dark halo log-slope to ±0.2 and normalization 
to ±15% at about twice the King radius. This is a factor 
of ~ 5 better than currently possible with LOS velocity 
dispersion data alone. The results from such observations 
will provide a very sensitive test of the CDM paradigm 
and an incisive tool for investigating the microscopic na- 
ture of dark matter. We estimate that with ~ 100 days 
of observing time over the 5 year lifetime of SIM, it will 
be possible to obt ain the proper motio ns of ~ 200 stars 
in multiple dSphs (jStrigari et al JI2007T ) . 
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